Non-Hermitian Hamiltonians with real eigenvalues coupled to electric fields: 
from the time-independent to the time dependent quantum mechanical formulation 
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We provide a reviewlike introduction into the quantum mechanical formalism related to non- 
Hermitian Hamiltonian systems with real eigenvalues. Starting with the time-independent frame- 
work we explain how to determine an appropriate domain of a non-Hermitian Hamiltonian and pay 
particular attention to the role played by PT-symmetry and pseudo-Hermiticity. We discuss the 
time-evolution of such systems having in particular the question in mind of how to couple consis- 
tently an electric field to pseudo-Hermitian Hamiltonians. We illustrate the general formalism with 
three explicit examples: i) the generalized Swanson Hamiltonians, which constitute non-Hermitian 
extensions of anharmonic oscillators, ii) the spiked harmonic oscillator, which exhibits explicit super- 
symmetry and iii) the — a; 4 -potential, which serves as a toy model for the quantum field theoretical 
4 -theory. 



I. INTRODUCTION 

In general most physicists will almost instinctively 
associate a non-Hermitian Hamiltonian with unstable 
states, decaying wave functions, resonances and dissipa- 
tion. Such type of systems have been studied for a long 
time. They arise for instance when coupling channels in 
a system in which the wavefunctions factorize into func- 
tions which depend on separate sets of variables. The 
effective Hamiltonians resulting in this manner are non- 
Hermitian and have complex eigenvalues [lj. However, 
one should note that Hermiticity of the Hamiltonian is 
only a sufficient condition, which guarantees real eigen- 
values and the conservation of probability densities. It 
needs to be emphasized that it is not a necessary con- 
dition and there could be non-Hermitian Hamiltonians 
with real discrete eigenvalue spectra, which then consti- 
tute potential candidates for physical applications, such 
as for instance atomic systems without decay. 

Precisely such type of Hamiltonian systems are cur- 
rently under intense investigation (for a collection of re- 
cent results see for instance [2|]). The central question 
in this context is of course how to obtain a consistent 
quantum mechanical framework. So far much effort has 
gone into the study of time-independent eigenvalue prob- 
lems. The main question we wish to address here is how 
to couple an external time-dependent electric field to a 
non-Hermitian Hamiltonian with real eigenvalues 

Our manuscript is organized as follows: For the bene- 
fit of the non-expert and the audience of this conference 
we commence in section II with a brief reviewlike intro- 
duction by recalling some by now well-known facts and 
arguments on the consistent quantum mechanical formu- 
lation of non-Hermitian Hamiltonian systems. In section 
A we review the time-independent formulation starting 
with a discussion of how to determine the appropriate do- 
main for a non-Hermitian Hamiltonian from the choice of 
asymptotic boundary condition. In part 2 of this section 
we explain the limited role played by 'PT-symmetry. In 
part 3 of section A we explain how pseudo-Hermiticity 
can be employed to map almost all relevant problems in 



the non-Hermitian scenario to a Hermitian system in the 
same equivalence class. The two systems obtained in this 
manner are therefore isospectral. In section B we discuss 
how this formalism can be extended to include an evo- 
lution in time. We describe here gauge transformations, 
perturbation theory and how to compute various phys- 
ical quantities in the non-Hermitian setting. In section 
III we discuss two methods of how to solve one of the key 
problems in this context, namely how to compute pseudo- 
Hermitian Hamiltonians. Section IV contains three ex- 
plicit examples to which the formulation from the previ- 
ous sections applies: i) the generalized Swanson Hamilto- 
nians, which constitute non-Hermitian extensions of an- 
harmonic oscillators, ii) the spiked harmonic oscillator, 
which exhibits explicit supersymmetry and iii) the — x 4 - 
potential, which serves as a toy model for the quantum 
field theoretical 4 -theory. We state our conclusions and 
an outlook to further problems in section V. 



II. THE GENERAL FRAMEWORK 

A. Time-independent quantum mechanical 
formulation 

1. The domain of non-Hermitian Hamiltonians 

The current interest in this subject was triggered eight 
years ago Q by the at the time rather surprising numer- 
ical observation that the Hamiltonian 



H = p — g(iz 



,N 



(1) 



defined on a suitable domain possesses a real positive and 
discrete eigenvalue spectrum for integers N > 2 with 
positive real coupling constant g. This property holds 
despite it being non-Hermitian H ^ H* and unbounded 
from below, for N — An with n G N . Throughout this 
paper we use atomic units % = e = m e — ca = 1. 

Viewing now H in |T]) as a differential operator in posi- 
tion space acting on some wavefunction &(z), one needs 
to specify appropriate boundary conditions in order to 
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select a meaningful domain. In [4] it was argued that 
the natural boundary condition, $(z) — > exponentially 
for \z\ — > oo, requires that one continues the eigenvalue 
problem into the complex z-plane. In fact, for H in ([1]) it 
was found that the wedges bounded by the Stokes lines 
in which this boundary condition holds are given by 



W L (N) = 
W R (N) = 



8 + N 
'2(N + 2) 
N 



7T < < 



N 



2(N + 2) 
4-N 



-TV < < 



(2) 



(3) 



2(N + 2) 2(N + 2) 

where 9 = argz. To see this one can follow the proce- 
dure for an asymptotic expansion of a differential oper- 
ator as outlined for instance in [5j]. Substituting &(z) = 
exp(<p(z)) into the eigenvalue equation iJ$ = e$ yields 
ip" + (ip r ) 2 + g(iz) N + e = 0. For \z\ — > oo with infinity 
being an irregular singular point one may assume that 
ip" <C {f 1 ) 2 ■ For large z we can also neglect e in compar- 
ison with the potential and obtain 

2^5 



ip(z) 



;(l+f) ? (l+f) 



AT- 



for 



(4) 



In order to extract the dominating exponential factor in 
<i>(z) and to achieve $(z) — » for \z\ — > oo, we require 
Re<p(z) < 0. With 6 = argz this is equivalent to 



TTiV 



N 



> 0, 



(5) 



which amounts to the conditions (|2|), (|3|) for the left and 
the right wedge Wl and Wr, respectively. Of course 
([5]) allows for many more solutions and therefore possi- 
ble wedges, but the selection criterion for @, Q is to 
reproduce the conventional wedge for the harmonic oscil- 
lator for N = 2, which is centered around the real axis. 

This means the domain of integration, which makes 
the eigenvalue problem of the non-Hcrmitian Hamilto- 
nian operator in (fT|) in position space well defined for the 
asymptotic boundary condition $(z) — * exponentially 
for |z| — > oo is any path in the complex z-plane which re- 
mains inside the wedges Wx and Wr when it approaches 
complex infinity. This means any path parameterized as 
z(x) with i£8, which satisfies 

(6) 



lim arg[z(x)] € W R/L 



guarantees the appropriate boundary condition, namely 
exponential decay at infinity of the wavefunction $(z). 
For various purposes, for instance when one is concerned 
about a fast numerical convergence, one can also deter- 
mine the anti-Stokes lines, that is the domain on which 
the wavefunction vanishes most rapidly, see e.g. Q. For 
H in (TTJ) the anti-Stokes angles 0^J R are just in the centre 
of W L and W R 0. 

Permissible domains are therefore usually some form of 
parameterizations for hyperbolae. For instance a modi- 
fied version of a parameterization used in Q was sug- 
gested in 0] 



zi(x) 



xcosl 



(7) 



with net This clearly satisfied the required asymptotic 
Jim^argMx)] = 9j*j L {N) S W R/L (N). (8) 

for all values of N. However, as we shall discuss in more 
detail in section IV certain manipulations depend cru- 
cially on the suitable choice of the parameterization and 
one needs various alternatives. The selection procedure 
for what is most "appropriate" is largely left to inspired 
guess work at this stage. As we shall see below, an ex- 
tremely useful variation of was provided in 



z 2 (x) = -2iVl + 



with 



lim arg[z 2 (x)] = -f G W R (N) 

x — >oo 

lim arg[z 2 (x)] = =f= G W L {N) 



(9) 



(10) 



for N — 3,4, ...,9. Another permissible parameteriza- 
tion can be found for instance [9j . We illustrate the above 
discussion with some examples in figure 1: 
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Figure 1: Stoke wedges in which the eigenfunctions of H in 
|T|) for iV = 4, 6 vanish exponentially when \z\ — ■> oo. Permis- 
sible paths z\ with a=l and Z2 as parameterized in and 
((9]), respectively. The Stokes lines are depicted as dotted lines 
in the figure. 
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2. VT -symmetry and real eigenvalues 

So how can one explain such unconventional behaviour 
that a non-Hermitian Hamiltonian possesses a real eigen- 
value spectrum? Shortly after the above mentioned ob- 
servation it was suggested that the reality of the spec- 
trum should be attributed to unbroken PT-symmetry 
ficj |. that is the validity of the two relations 

[H,VT]=0 and PT$ = $, (11) 

where $ is a square integrable eigenfunction on some 
domain of H. In other words, when the Hamiltonian and 
the wavefunction remain invariant under a simultaneous 
parity transformation P and time reversal T 

V : p — > — p z — > — z 

T : p — > — p z — > z i — > —i (12) 

VT : p — > p 2 — > — z i — > — i, 

the eigenvalues of H are real. As an example one sees 
that obviously the Hamiltonian in equation {lj is PT- 
symmetric. What is less straightforward to see is that 
for N < 2 the second relation in (fTTj) does not hold. 
Analytic arguments, which establish these facts for the 
Hamiltonian fTJ may be found in [n], [12] ■ 

We shall now outline to what extend PT-symmetry 
can be utilized. Clearly V 2 = T 2 = (VT) 2 = I and 
the last relation in (JTSJ) implies that the PT-operator is 
an anti-linear operator, i.e. it acts as VT(X^> + n^>) = 
A*PT$ + p*VT^> with A, /i e C and * being eigen- 
f unctions of the Hamiltonian H with eigenenergies e, 
HQ = £$. The anti-linear nature of the PT-operator 
serves well to establish the reality of the spectrum, 
i.e. e = e*, when both relations in pip hold. This follows 
simply from 

e$ = H$ = HVT<$> = VTH<£ = VTe<5> = e*VT<£ 
= (13) 

Unfortunately, the anti-linearity is also responsible for 
the possibility that only the first identity in (|11[) could 
hold, but not the second. In this situation one speaks of 
a broken PT-symmetry. The argument leading to this is 
straightforward [HI, [3| : Let us consider first a unitary 
operator U for which by definition 

(Z7*|f/$) = (#|$) (14) 

holds for all eigenfunctions $, of H . From equation 
(ITU) it follows that = with |u| = 1 for all which 
means that a unitary operator has only one dimensional 
representations. This property changes for anti-unitary 
operators A, as in that case only A 2 is a unitary operator, 
which can be seen from 

(A 2 ^\A 2 Q) = (A$\A$!) = (#|$) . (15) 

Now we can only deduce from (|15p that A 2 ^ = a 2 ^ with 
a 2 = 1 for all and this means that an anti-unitary 



(which is implied by anti-linearity) operator could have 
a two dimensional representation A^ — a*$, A$> = a^. 
Indeed when a is purely imaginary one can not construct 
a linear combination £1 = A$ + p^>, with A, /x € C of the 
two so-called flipping states <&, ^P, which remains invari- 
ant under the action of A. We see that 

Afi = A*a* + pV$ = A$ + ^ (16) 

implies that /i = A*a, A = p*a* and therefore a 2 = 1. 
This means that only for a = ±1 the two-dimensional 
representation is reducible and for purely complex a it 
is irreducible. In the latter situation the second relation 
in (fTTj) does therefore not hold. From (|13p we see that 
PT$ is an eigenfunction of H with eigenvalue e* when 
$ is an eigenfunction of H with eigenvalue e. Thus when 
the second relation in (fTTj) does not hold the eigenvalues 
of H come in complex conjugate pairs. 

Thus PT-symmetry is merely a fairly good guiding 
principle and serves to identify immediately potentially 
interesting non-Hermitian Hamiltonian systems. How- 
ever, as argued above the PT-symmetry of H does not 
constitute a guarantee for a real eigenvalue spectrum. 
It remains an open question at this stage to determine 
under which circumstances the PT-symmetry is broken, 
albeit for Hamiltonians acting in a finite dimensional 
Hilbert space an algorithm based on stability theory has 
been provided [l5|. In addition one should stress, that 
PT-symmetry can not be regarded as the fundamental 
property, which explains always the reality of the spec- 
trum for non-Hermitian Hamiltonian systems as there ex- 
ist also examples with real spectra for which not even the 
Hamiltonian is PT-symmetric 0, [H| (see also examples 
below). In fact, more fundamental is the necessary and 
sufficient condition that the Hamiltonian must be Hermi- 
tian with regard to some positive definite inner product 
[13] as we shall discuss next. 

We summarize the role played by PT-symmetry in fig- 
ure 2: 




Figure 2: PT-symmetry and real eigenvalues 
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3. Pseudo-Hermiticity and real eigenvalues 

The formal question of how to establish a consis- 
tent quantum mechanical formalism for non-Hermitian 
Hamiltonian systems has already been discussed in [l8| 
prior to the above mentioned numerical observation. In 
fact the possibility of extending a Hilbert space by a new 
intermediate state, which then leads to an indefinite met- 
ric has already preoccupied particle physicists more than 
half a century ago 19]. Some of these old results have 
recently been re-discovered and developed further. As 
already mentioned Hermiticity is a useful property as it 
guarantees the reality of the spectrum. Let us briefly re- 
call this standard argument and discuss how it needs to 
be altered for the present scenario. 

Suppose we have a diagonalizable Hermitian (symmet- 
ric) operator h with regard to the conventional inner 
product 



{<t>n \h(f>m) = (h4>n\(l>r, 



(17) 



We use here Hermiticity in the sense that it implies self- 
adjointncss and ignore possible subtleties, which might 
arise from domain issues. In general we understand here 
the domain to be the entire real axis. Multiplying next 
the eigenvalue equations 



\h(j) m ) = e m \<f> m ) 



and 



{h<t> n \ = e* n {ct> n \. (18) 



by (<p n \ and \(f> m ), respectively, we obtain 

{4>n \h(f>m) = £m {4>n \4>m) 
{h(f> n \4> m ) = £* {(j) n \4>m) ■ 



(19) 
(20) 



Taking the difference between p^|) and (|2"0)) thus implies 
for n = m that Hermiticity of h with regard to the stan- 
dard positive-definite inner product (4> n \4>m), i-e. the va- 
lidity of (fl7|) , is a sufficient condition for the energies e n 
to be real. Taking next n ^ m Hermiticity then also 
implies the orthogonality of the states \<p n ) for all n. 

It turns out that for non-Hermitian operators we 
only need to change the definition for the inner prod- 
uct, i.e. change the metric, to draw the same conclu- 
sions [13, EH- Taking now the domains as discussed 
in section 3, by definition we obviously no longer have 
(§ n \H§ m ) = (H$ n \& m ) for a non-Hermitian operator 
H with <E>„ obeying the eigenvalue equation 



|4>, 



(21) 



Therefore there is no guarantee for the reality of the spec- 
trum and neither for the orthogonality. However, assum- 
ing rj to be a Hermitian operator with respect to the 
standard inner product, we can define a new inner prod- 
uct 



($n|?/ 2 $m) ■ 



(22) 



Supposing now that H is Hermitian with regard to this 
new inner product 



we may employ exactly the same arguments as above 
and ensure the reality of the spectrum as well as the 
orthogonality (Qnl&m}^ = <5„ im . Note that with regard 
to the standard inner product one finds in general that 
($n |$m) + &n,m, see e.g. [ltj. 

What is left is to characterize in more detail and possi- 
bly to determine is the metric operator rj 2 . Mostafazadeh 
US HH HH proposed to assume that H is a pseudo- 
Hermitian operator satisfying 

h = ijHtj- 1 = h ] = rf 1 ^!) H*=r?Hrr*, (24) 

where rj is the Hermitian operator with regard to the 
standard inner product as introduced above. Since the 
Hermitian Hamiltonian h and the non-Hermitian Hamil- 
tonian H are related by a similarity transformation, they 
belong to the same similarity class and therefore have the 
same eigenvalues. The corresponding time-independent 
Schrodinger equations are then simply (fT5)) and (|2ip . 
where the corresponding wavefunctions are related as 



$ = rj <f>. 



(25) 



Having real eigenvalues for the Hermitian Hamiltonian 
h then guarantees by construction the same real eigen- 
spectrum also for H. In fact the necessary and sufficient 
condition (|23p , which ensures the reality of the spectrum 
for H follows then from (fTT|) 



•q 1 tp\rj 2 Hrj (; 



(#|-H"$) = {^\rj 2 H<S>) = 
= (iPlriHr/- 1 4>) = {ip\h(j>) = {hip\<f>) = {<qHr)- l tl)\(f>) 
= (H^\t]4>) = (HV\ ?y 2 $) = (HV\ $), ( . (26) 

Clearly, when -q is Hermitian with regard to the stan- 
dard inner product it is also Hermitian with regard to 
the r/-inner product (|2"2"]) . 

A particular example for an ?7-inner product is the one 
introduced in [lOl 



m) CVT := (CTT\*)) 7 



|*) , (27) 

where the new operator C, with C 2 = I, [H,C] = and 
[C,PT] = 0, is employed. In position space it reads 
C(x,y) = J2n^ n ( x )^ n (v)- The operators C and rj 2 are 
simply related as C = r\~ 2 V . 

In addition one should stress that in fact these inner 
products can also be derived [H, [24[ when starting from 
a biorthonormal basis, which is quite common to use in 
the study of non-Hermitian Hamiltonian systems with 
complex eigenvalues, that is decaying states, see e.g. [25|, 

m. 

Crucial for a proper quantum mechanical framework 
is of course to clarify the nature of the physical observ- 
ablcs. In order to be suitable for a physical interpreta- 
tion, observables O have to be Hermitian operators act- 
ing in some physical Hilbert space. From what has been 
outlined above with regard to the inner products, it is 
natural to take them to be Hermitian with respect to the 
new 77-inner product 



(23) 



(28) 
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This implies immediately that when o is an observable in 
the Hermitian system, then 

O = r)- l aq & O = r/^O^r] 2 (29) 

is an observable in the non-Hcrmitian system. This 
means in turn that the standard position operator x and 
the momentum operators p are in general not observ- 
able in the non-Hcrmitian system, but rather their non- 
Hermitian counterparts X and P, respectively. Clearly 
X and P satisfy the standard canonical commutation re- 
lations [X, P] = i when [x,p] = i. For Hamiltonians h, H, 
which admit a polynomial expansion in {x,p}, {X, P}, 
it follows then directly from (|24[) that 

H{x,p) = ri- l h{x,p)r) = h(X, P), (30) 
h(x,p) = r)- 1 H' f {x,p)<ri = H Jf {X,P). (31) 

These relations serve for instance as a consistency check 
when we start with a given non-Hermitian Hamiltonian 
and construct its Hermitian counterpart by means of a 
similarity transformation. Moreover (|30|) provides a sim- 
ple way to express the non-Hcrmitian Hamiltonian in 
terms of the canonical (X, P)-variables, which have a 
physical meaning for that system rather than the (x,p)- 
variables, which arc in general meaningless in that con- 
text. In addition, one may use (|3"0|) as a principle to 
construct non-Hcrmitian Hamiltonians with real spectra 
from a given Hermitian Hamiltonian and a set of canon- 
ical variables and vice versa with (|3ip . When rj is VT- 
symmetric in the (ir,p)-variables the corresponding quan- 
tities in the non-Hermitian system will be PT-symmetric 
in the (X, P)-variables. 

We conclude with a final comment in regard to the 
uniqueness of the metric operator rj 2 . In fact there are 
various types of ambiguities arising, which we comment 
on in section III when we discuss how to compute rj 2 ex- 
plicitly. In [l8| Scholtz, Geyer and Hahne proved that the 
metric operator rj 2 is uniquely determined on a Hilbert 
space if and only if a set of observables Oi with respect to 
(|28p is irreducible on this Hilbert space. The latter means 
that there is no bounded operator besides the identity, 
which commutes with all observables Oi. Taking this 
result into account allows to move the nature of the am- 
biguities from the metric to the specification of the set of 
observables. As we shall see below a subset or even one 
observable might be enough in practise. 

B. Time-dependent quantum mechanical 
formulation 

Let us now discuss how to couple a laser field to the 
non-Hcrmitian Hamiltonians H , which have the proper- 
ties described above. In the simplest scenario, i.e. if the 
parameters involved lie within a non-relativistic regime 
and the dipole approximation holds, such a field can be 
approximated by a time-dependent electric field E(t). In 
the following, we will briefly recall our recent results on 



the temporal evolution of the resulting system [3j . For 
simplicity, we assume that E(t) is linearly polarized and 
has a finite duration r. 

Within a Hermitian framework and in the length 
gauge, such an evolution is described by the time- 
dependent Schrodinger equation 

id t <f>{t) = hi(t)</>(t), (32) 

where 

2 

hi (*) = y + V(x) + xE(t) =h + xE(t) (33) 

is the Stark-LoSurdo Hamiltonian [27|, HU . For a pulse 
of finite duration, h(j)(0) = £</>(0) and /i0(r) = £</>(r). Wc 
assume here that h possesses a non-Hermitian counter- 
part H which is in the same equivalence class, i.e. the 
validity of the first relation in (|2"4")l . 

1. Time-evolution operators 

The central quantity of interest in this context is the 
time-evolution operator 

u(t,t') = Texp dsh(s) \ , (34) 

which evolves a wavefunction from a time t' to t, that 
is <j){t) = u(t,t')(f>(t'). In |[3"I|). T denotes the time or- 
dering. One should note that, in general, u{t, t') ^ 
exp [— ih(t — t')] . In fact, such a relation only holds for 
Hamiltonians which are not explicitly time-dependent, as 
is not the case for the scenario we have in mind. 

When h(s) is a self-adjoint operator in some Hilbert 
space, u(t,t') satisfies the relations (29l. l30l. [3TI] 

id t u(t,t') = h(t)u(t,t'), 
u(t, t')u(t', t") = u(t, t") and u(t, t) =1 . (35) 

We will now assume that the similarity transformation 
rj extends to the time-dependent case. Thus, H(t) — 
r}~ x h(t)r] , with H(t) ^ H\t). We take rj to be time- 
independent. This allows us to guarantee that the rela- 
tions 

id t U(t,t') = H(t)U(t,t'), 
U(t, t')U(t', t") = U(t, t") and U(t, t) = I, (36) 

for the time-evolution operator U(t,t') = rj~ 1 u(t,t')r] as- 
sociated to the non-Hcrmitian Hamiltonian H(t), also 
hold. Then this operator fulfills the condition W(t, t') = 
rfU- 1 ^,^- 2 , which follows from u+ (*,*') = u _1 (t,i / ). 

It is worth stressing that we make no simplifying as- 
sumption on the time dependence of the Hamiltonian. In 
fact, the only requirements involved in our approach are 
i) that H be pseudo-Hermitian and ii) that the similarity 
transformation rj be time- independent. Such conditions 
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guarantee that the time-dependent Schrodinger equation 
and the relations involving the time evolution operator 
remain valid also in the non-Hermitian case. 

These conditions, however, are far more general than 
those normally encountered in the literature. In fact, 
most studies make several simplifying assumptions on 
the time-evolution operator, in the sense that they con- 
centrate on Hamiltonians which are either not explic- 
itly time-dependent, or which vary adiabatically and/or 
periodically with time. The first scenario is addressed 
by either solving the eigenvalue problem H<& = e$, 
or, at most, by employing the time-evolution operators 
U(t,t')=exp[-iH{t-t')]. 

The remaining situations are widespread in the atomic 
physics literature, in the context of open quantum sys- 
tems. Roughly speaking, if a system is close to the adia- 
batic limit this means that it is varying so slowly that the 
problem can be reduced to solving eigenvalue equations of 
the form H(t)Q(t) = £„(t)$(t). In a standard, Hermitian 
framework, this implies that dtu(t)u> (t) <C u(t)h(t)u'(t), 
and that transitions between different time-dependent 
eigenstates of H(t) will be induced by perturbations 
around the adiabatic limit. Specifically for a system 
coupled to an external laser field, the time-dependent 
energies e n (t) give the field-dressed states (for a first 
derivation of the adiabatic theorem and for an extension 
of such a theorem to non-Hermitian open quantum sys- 
tems, see dH and [Hj|, respectively). For periodic fields, 
such a procedure is closely related to the Floquet theory, 
for which there also exists time-dependent "quasiener- 
gies". This approach may be problematic if the field 
varies abruptly with time, such as, for instance, if it is 
an ultrashort pulse. 



2. Time-dependent physical quantities 

The time-evolution operators characterized in the pre- 
vious subsection may then be employed to compute var- 
ious quantities of physical interest, such as for instance 
the transition probability 



$„| U(t,0) $ m ) = \(4> n \u(t,0) cj> m )\ 2 , (37) 



from an eigenstate \4> m ) to \(f> n ) of the Hermitian elec- 
tric field-free Hamiltonian h or eigenstate |$ m ) to 
of the non-Hermitian electric field-free Hamiltonian H . 
Another physical quantity of interest is the time evolu- 
tion for the expectation value of an observable in the 
state n is 



O n (t) = (£/(i,O)$„(O)|0f/(t,O)$„(O)) 7 

= (u(t,0)&,(0)|ou(i,0)&.(0))„ 
= o n (t). 



(38) 



the relevant computations for the non-Hermitian system 
can be translated into the Hermitian one, provided the 
77-operator is known. 



3. Gauge transformations 

Apart from employing the length-gauge Hamiltonian 
hi(t), one may describe a Hermitian Hamiltonian system 
coupled to an electric field in other gauges. Concrete 
examples are the velocity gauge, obtained by employing 
the minimal-coupling prescription p — > p — b(t), or the 
Kramers- Henneberger gauge, obtained with the shift x — > 
x — c{t) in the field-free Hamiltonian h as introduced in 
(|33p . The corresponding Hamiltonians are given by 

hv(t) = (P ~y ))2 + V(x) = h(p - b(t)) (39) 



and 



h KH (t) = ^- + V(x- c(t)) = h(x - c(t)), (40) 



respectively. In equation (|3"9"|) and ([4*U|) . 



b(t) = / dsE(s), c(t) = / dsb{s) 



In a similar way we may proceed to compute ionization 
rates and probabilities etc., but these examples are suf- 
ficient to see that, as in the time-independent scenario, 



(41) 



are the momentum transfer b(t) from the laser field to the 
system in question and the classical displacement c(t) in 
the system caused by the laser field. 

Depending on the problem at hand, the gauge choice 
may considerably facilitate the computations. For in- 
stance, the length gauge is very appropriate for pertur- 
bation theory in the electric field, as the field coupling 
involves only one additional term, or for physical inter- 
pretations in the low-frequency regime, since it allows the 
physical picture of an effective time-dependent potential. 
The Kramers-Henneberger gauge is most useful in the 
high-frequency regime, especially if one wishes to exploit 
the periodicity of the field and perform Floquet expan- 
sions. Each formulation can be obtained from the other 
employing gauge transformations. The Hamiltonians in 
the length, velocity and Kramers-Henneberger gauge are 
related by 

hi(p,x)-xE{t) = h v (p+b{t),x) = h KH (p,x+c{t)). (42) 

We will now perform such transformations for non- 
Hermitian Hamiltonian systems. First, we will replace 
the wavefunction cf> in the time-dependent Schrodinger 
equation related to the Hamiltonian h by <p = a{i)~ l 4>' , 
with a(t) being some unitary operator. This yields 
[Il,[IJl3l| 

ift0' = V(t)^'=[o(t)/i(t)o(t)- 1 +ifto(t)o(t)- 1 ] (j)'. (43) 

Specifically, the standard transformation from the length 
to the velocity gauge, and from the velocity to the 
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Kramers-Henneberger gauge, which are extensively used 
in strong-field laser physics, are given by, 



at. 



„(t) = e lb{t)x and a v ^ K ii{t) = e^We^. (44) 



respectively. In equation (I44p , in addition to the momen- 
tum transfer and classical displacement, we have also in- 
troduced the classical energy transfer d(t) — | f Q dsb(s) 2 . 

If the system is pseudo-Hcrmitian, one may employ the 
relation <p = 77$ to obtain the gauge transformation 

i d t $'=[A(i)iJ(t)A(i)- 1 +id t A{t)A{t)- 1 } (45) 

where 

a(t) = rjAit)^ 1 and h(t) = •qH(t)r)~ l , (46) 

and the expression in brackets, on the right-hand-side of 
(|4"B")) . denotes the gauge-transformed Hamiltonian H'(t). 
The gauge transformations A(t), as it should be, guar- 
antee the invariance of the physical observables, when 
computed using the generalized inner product (|22p . Now 
the relations 

H l {X,P) = H v {X,P + b{t))+XE{t) (47) 
= H KH {X + c(t),P)+XE(t), 

hold for pseudo-Hermitian Hamiltonians. 

4- Perturbation theory 

Since, in most realistic situations, the time-dependent 
Schrodinger equation describing the evolution of a sys- 
tem with a binding potential V(x) subjected to a time- 
dependent laser field E(t) does not possess an analytic 
solution, it is necessary to resort to perturbation theory. 
In order to construct a perturbative series in a pseudo- 
Hermitian framework, we will initially consider a time- 
dependent Hermitian Hamiltonian h(t) = ho(t) + h p (t), 
where ho(t) and h p (t) are also Hermitian and satisfy 
the time-dependent Schrodinger equation. Using the Du 
Hamel formula [29|, [3(| HH, we can express the time- 
evolution operator u(t,t') associated to h(t) as 



u(t,t')=u (t,t') 



1 J u(t, s)h p (s)uo(s,t')ds, 



(48) 



where Uo(t,t') is the time evolution operator with re- 
spect to ho(t). Equation (|4*8"|) can then be solved itera- 
tively to an arbitrary order in h p (t), which will be the 
perturbation. Roughly speaking if h p {t) <C ho(t), the 
series obtained by such means has a great chance to con- 
verge. For instance, for weak laser fields and in the length 
gauge, a natural choice is to take h p (t) = xE{t) and 
h (t) = p 2 /2 + V, whereas in the strong-field regime we 
take ho(t) — p 2 /2 + xE{t) as the Gordon- Volkov Hamil- 
tonian and the perturbation is chosen as h p (t) = V. 

Similarly, for the time evolution operator U(t,t') 
related to its pseudo-Hermitian counterpart H(t) — 



H (t) + H p (t), with H (t) = J7~%(i)»7 and H p (t) = 
■q h p {t)rj, we may also write 

U(t, t') = U (t, t')-i ( U(t, s)H p (s)U (s, t')ds, (49) 



where Uo(t,t') is related to the Hamiltonian Ho(t), and 
the perturbative series is obtained by iterating equation 
(|4*^| up to the desired order. 

a. The weak intensity regime As argued in the pre- 
vious subsection one can in general not compute the time- 
evolution operator exactly and has to resort to perturba- 
tion theory instead. We illustrate here briefly how this 
works more explicitly in the different intensity regimes. 

We commence with the weak intensity regime and we 
will consider first-order perturbation theory with respect 
to the external laser field amplitude E . Iterating ([4*5]) 
it follows that to this order the time-evolution operator 
can be approximated by 



U (1) (t,0) = U Q (t,0)-i / U o (t,s)XE(s)U Q (s,0)ds, (50) 



where Uo(t,0) — exp[—iHt]. Subsequently the transition 
probability (|37[) from a state m to n to this order becomes 



T~*n< — m — 



(*n|*$ f . 



dse l{e 



'2500 



(51) 



Note here the occurrence of the matrix element 
(^ n \X^ m } v — (</>„| xcf) m ) , which results from taking the 



non-Hermitian version of the Stark-LoSurdo Hamiltonian 
in ([33J) to be Hi{t) = H + XE{t). In case we add xE(t) 
instead of XE(t) the amplitude {<j) n \ nxrj^ 1 <p m ) would 
occur. With our examples below we demonstrate that 
the latter matrix element is very often unphysical. 

b. The strong field regime Next we will address the 
opposite scenario, namely the situation in which the laser 
field is larger, or at least comparable to the atomic bind- 
ing forces. Such a physical framework has become of in- 
terest since the mid-1980's, when intense lasers became 
feasible, due to the wide range of phenomena and ap- 
plications existing in this context. Concrete examples 
are high-order harmonic generation, above-threshold ion- 
ization, or laser-induced single and multiple ionization 
(for reviews we refer to [13, HU, HI[). In this case, it 
is a common procedure to perturb around the Gordon- 
Volkov Hamiltonian, which, in a non-Hermitian frame- 
work and in the length gauge, is given by h[ GV ^ (t) = 
P 2 /2 + XE(t). To first order, the time-evolution opera- 
tor then reads 



uW(t,o) = u{ GV \t,o) 



(52) 



i / U^ GV \t,s)V(X)U^ GV \s,0)ds, 
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where the Gordon- Volkov time-evolution operator is 
given by 

u[ GV \t,0) = A^/Wcxp[-zP 2 i/2]A A -^ ; (0). (53) 

The gauge transformation AKH^i{t), from the Kramers 
Hcnneberger to the length gauge, is written as 



with ho = ho, hi = h\. For the case when one has 
the condition c^ +1 \ho) — for some finite integer t, we 
found in [3j the closed expressions 



n=l v ; 



(59) 



Akh^i (t) = r 1 - 1 e lc(t)p e ld ^e' lb{t)x r 1 (54) 

(55) 



= & ic{t)P e id(t) e -ib(t)X _ 



Obviously, one may also define a Gordon- Volkov Hamil- 
tonian in the velocity gauge as H[ GV \t) — (P — b(t)) 2 /2. 
In this case, the corresponding time evolution operator 
is t4 Gy) (i,0) = e tb ^ x ul; GV) {t,0)e-^ x . 



III. COMPUTING PSEUDO-HERMITIAN 
HAMILTONIANS 

Having discussed the central role played by pseudo- 
Hermitian Hamiltonians it is vital to have a constructive 
method to realize them. In other words we wish to com- 
pute Hamiltonians h = h' and H ^ H> belonging to the 
same equivalence class. This is a well defined problem, 
but in most cases very difficult to solve. Here we present 
two different types of methods to achieve this. 



Similarity transformations from operator 
identities 



Supposing that the similarity transformation (|24p can 
be realized using a Hermitian operator of the form i] = 
exp(g/2), the second relation in (124| implies by standard 
Baker-Campbell-Hausdorff commutation relations that 

Hi = H+[q,H} + ± [q, [q, H}} + i [q, [q, [q, H}}} + ... 



3! 



(56) 



ra=0 



For convenience we have introduced here a more compact 
notation for the n-fold commutator of the operator q with 
some arbitrary operator O as 



4")(0) := [q, [q, [q, . . . [q, O] . . 



(57) 



Clearly, if for some integer n the n-fold commutator 
Cq^ (H) vanishes the conjugation and therefore the sim- 
ilarity transformation can be computed exactly. In or- 
der to see this more explicitly we separate next the non- 
Hermitian Hamiltonian into its real and imaginary part 
and bring it into the form 



ff^o-E^^cf^CM, (60) 

which are related according to the first identity in (I24|) . 
Here [x] denotes the integer part of a number x. The E n 
are Euler's numbers 

E x = 1, E 2 = 5, E 3 = 61, E 4 = 1385, . . . (61) 

and the K2n-x may be computed from them according to 



1 I(»+i)/2] 



/•■„ = - > (-l) n+m ( n \E„ 



2" ^ 

m— 1 



\2m 



The first examples arc 



«i = 2' K 3 



1 1 



17 



(62) 



(63) 



Depending on how large I becomes the explicit evalua- 
tion of sums in (|59p and (|60|) can become rather compli- 
cated. In fact, in most cases the series does not terminate 
and one has to compute the expressions perturbatively. 
We shall not discuss such cases here and refer instead to 
the literature [3, H3, [H, [H, 5(| ■ 



B. Similarity transformations from differential 
equations 

Alternatively one can follow a proposal put forward by 
Scholtz and Geyer jJJ, |42j and solve ([24]) by means of 
Moyal products instead of computing commutators. The 
central idea is to exploit isomorphic relations between 
commutator relations and real valued functions multi- 
plied by Moyal products, which correspond to differen- 
tial equations. We shall demonstrate that this approach 
is rather practical and allows to compute pairs of isospec- 
tral Hamiltonians h — h) and H ^ H\ when they are of 
polynomial nature. 

We use a slightly different definition for the Moyal 
product as in 41. 42], since then the resulting differen- 
tial equations become simpler [3{| . Following for instance 
[IH we define the Moyal product of real valued functions 
depending on the variables x and p as 

f(x,p)*g(x,p) = f(x,p)ei<tt-tt)g(x,p) = (64) 



H = ho + ih±, 



(58) 



EHrEt- 1 )*! * )9 t x d;- t f(x,p)d s x -%g(x,p). 
s =o ' t=o 
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One may then use this expression to turn all operator 
identities into differential equations. In principle this 
yields differential equations of infinite order, but when 
f(x,p),g(x,p) are of polynomial nature the series termi- 
nates and the order will be finite. For instance, if we want 
to compute the commutator [x,p] — i we have to evaluate 
the corresponding Moyal product relation x*p~p*x — i. 
Here and in some places below we emphasize the oper- 
ator nature of the quantities involved by dressing them 
with hats. In order to keep notations simple we do not al- 
ways make this rigorous distinction, when it is not strictly 
necessary. Matters become more complicated when the 
resulting real valued function depends on x as well as on 
p. As for a function the ordering is of course irrelevant 
we need a prescription of how to turn such a function 
back into operator valued expressions. Computing for 
instance 



Aipx 



2 2 

x *p 



2 2 
p * X 



Aipx, (65) 



we observe that we obtain the correct operator valued 
expression for the last equality when we replace px — > 
(px+xp)/2. In general we have to replace each monomial 
p m x n or x n p m by the totally symmetric polynomial S m ,n 
in the m operators p and n operators x 

'" ! "' ^ (66) 



S m ,n , , / 

[m + ny. ^— ' 



^p m x n . 



The sum extends over the entire permutation group ir. 
For our purposes we have usually a given non-Hermitian 
Hamiltonian H and wish to compute from the second 
relation in (|24| the Hermitian operator rj 2 . The corre- 
sponding differential equation is then simply 



H^(x,p) -kr] 2 (x,p) = i] 2 (x,p)-kH(x,p). 



(67) 



Subsequently, one may compute also rj(x,p) and h(x,p) 
in a similar manner. 

A comment is due concerning the uniqueness of the so- 
lutions. Having solved various differential equations, we 
naturally expect some ambiguities in the general solu- 
tions, which mirror the possibility of different boundary 
conditions. However, one should emphasize that these 
ambiguities are not only present when using Moyal prod- 
ucts, but are a general feature occurring also when using 
commutation relations of the type (|59|) and (f60|) . It is 
clear that in that context one may only fix the operator 
q up to any operator which commutes with the Hermi- 
tian part of H, that is ho- This means that, in (|59|) and 
(I60p . the expressions are insensitive to any replacement 
q — > q + q with [q, ho] = 0. A further type of ambigu- 
ity, which is always present is a multiplication of rj 1 by 
operators which commute with H , i.e. we could re-define 
i] 2 — > rj 2 Q for any Q, which satisfies [Q, H] = 0. 

It should be mentioned that there are also other pos- 
sibilities to evaluate the similarity transformations, such 
as for instance suggested in [44| or directly by using prop- 
erties of differential equations [45| . 

Let us now demonstrate with some concrete examples 
how the above mentioned formalism can be applied. 



IV. (QUASI) EXACTLY SOLVABLE MODELS 

Non-Hermitian Hamiltonians may arise for various dif- 
ferent reasons. In the following we provide three such 
examples, which all arise from quite different argumen- 
tations and thus provide several types of motivations to 
study non-Hermitian Hamiltonian systems. 



A. The generalized Swanson Hamiltonian 

One type of non-Hermitian Hamiltonian system arises 
form a purely mathematical consideration simply by 
perturbing a Hermitian Hamiltonian by adding a non- 
Hermitian term. We start with a straightforward ex- 
ample, which results when perturbing the anharmonic 
oscillators 



h° n (a) = l P 2 + ~x« 



(68) 



for n = 1,2,3,... and a£i Defining now the Hermitian 
operators rj m = exp(q m /2) with q m = 2g/mx m for m — 
1,2,3,... it is straightforward to compute that 

c£ (/£(«)) = ig(px m - x + x m - l p) (69) 
c£>(A°(a)) = -4 5 V™- 2 (70) 
cfj«(a)) = (71) 

for all n, m > 0. With (f6"9"j) - (j7"Tj) the generic expressions 
O and (JBDJ) yield with I = 2 



hC s m (a,g) = h° n (a) + ~g 2 x 2m - 2 



(72) 



H^i(a,g) = h° n (a) - i^px™' 1 + x^p), (73) 



which are related according to the first relation in 
In the special case n = m = 2, the Hamiltonian Hf? § be- 
comes the Swanson Hamiltonian discussed in [42|, Ey, [47j 
upon some change in the conventions for the coupling 
constants. This Hamiltonian arises in the second quanti- 
zation H = ciaa+C2a' ! +C3d^ a where the Cj are coupling 
constants and a) = {x — ip)/y/2, a = (x + ip) /v2 are the 
usual creation and annihilation operators, respectively. 
The sequence of Hamiltonians (|75)l illustrates our asser- 
tion on the limitations of 'PT-symmetry in section II A 1, 
that there are non-Hermitian Hamiltonians with real en- 
ergy spectra which are, however, not 'PT-symmetric. As 
one easily sees H'^f n {a 1 g) is not 'PT-symmetric when m 
is odd, but still has a Hermitian counterpart and there- 
fore real eigenvalues. 

Let us next assume that we had simply given the non- 
Hermitian Hamiltonian and we wanted to compute the 
r/-operator. For instance, for i/^f (a, g) and -ff^pf (a, g) 
the corresponding equations (f6"T)) become 



= Agpxif + 2axd p T] 2 - 2pd x r] 2 + gd p d x rj 2 , 
= igpxi] + 4:ax 3 d p r] 2 — 2pd x rj 2 + gd p d x rj 2 — axd^rj 



(74) 

3„2 
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respectively. Both equations are easily solved by rj 2 = 
exp(gx 2 ), thus confirming our previous calculation. 

Having the operator rj — exp(gx m /m) at hand we com- 
pute from (|29p the observables which correspond to the 
position and momentum operator in the non-Hermitian 
systems H^f n {a,g) as 

X = x and P = p-igx m ~ 1 , (75) 

respectively. Then it is easily verified that indeed ([50]) . 
(|31|) and g7D hold. 

With regard to the uniqueness of this solution one 
can see that the first equation in (f71| is also solved 
by fj 2 — exp(—g/ap 2 ). In fact for what has been re- 
marked at the end of the last section, it is clear that 
there should be more solutions corresponding to ff 2 — 
exp(gx 2 ) / (h^iaj) , with / being some arbitrary well be- 
haved function restricted by the demand that $ = 77" 1 
remains a bounded function. Obviously fj 2 — fj 2 for 
f(x) — exp(— 2gx/a). To see that other choices for f(x) 
will also lead to solutions of ([7i]) is less straightforward 
as we have to turn the operator valued expressions for ff 2 
first into real valued functions before we can verify (|74|) . 

Let us next illustrate how to fix the ambiguities by an 
explicit choice of the observables in the non-Hermitian 
system, which is always possible for what has been said 
at the end of section II A. Demanding for instance that 
X = x should be an observable in the non-Hermitian 
system, it follows immediately that the only choice for 
f{x) is f(x) = 1 and therefore (|75|) is the corresponding 
set of canonical variables. In turn we could also choose 
P = p to be an observable, which leads to fj 2 and X = 
x — ig/ap. For m 7^ n it is not possible choose P = p 
to be an observable as one can not find a function f(x) 
such that fj 2 becomes a function of p only. 



B. The spiked Harmonic Oscillator 

A further interesting example is the spiked harmonic 
oscillator as it exhibits an explicit supersymmetry (48l . 
|49L [50j and therefore also phenomena like degeneracy of 
the energy eigenvalues and even level crossings. The Her- 
mitian version of this Hamiltonian is simply 

h SHO {x,p) = \ V 2 + X 2 x 2 + "Lllll. (76) 
2 ar 

This example is very instructive as it is exactly solvable. 
The normalized eigenfunctions are 

= (-ir M^^ e-^x^L^Xx 2 ), (77) 
y 1 (a + n + 1) 

where the L"(x) denote the generalized Laguerre poly- 
nomials and the eigenenergies are 

el = A(4n + 2a + 2). (78) 



Clearly there is a degeneracy of the energy levels for 
e~ a = £n- a - The standard harmonic oscillator Hamilto- 
nian results from ([76| for a — ±1/2. The corresponding 
wavefunctions are related to (|77|) as tfil/ 2 = * 2n— 1 02n+li 

(jhi 1 ^ 2 = (— ^) n( t ) 2n > ■ The motivation here to introduce 
an Hermitian counterpart for this Hamiltonian is that 
one wishes to regularize the singularity of the potential 
at x = 0, see e.g. [50j | . 

With 77 = exp(— £p) one easily produces the desired 
shift and with (|24|) one obtains 

H SH °(x,p) = \p 2 + X 2 (x - zO 2 + f- (79) 
2 (x- iQ 2 

This is an example for which the Moyal products are not 
very suitable for the computations as the last term in 
the potential of (I79p is responsible for the fact that the 
related differential equations are of infinite order. 

Nonetheless, commutators are easily evaluated in this 
case and for instance the canonical variables for the non- 
Hermitian system are computed in a rather trivial way, 
resulting to 

X = x - i£ and P = p. (80) 
Once again we verify ([30]) and (|3Tj) for consistency. 



P(co) 







1 1 1 1 


1 ' 




x = 20 nl to 


,■ \ 






5=0 

5 = 2 

5 = 5 

1= 10 








t 

1 

1 

\ 

1 

1 








,' 1 

1 ; 
: 1 
1 *, 

1 ' 

• .' ■ 1 
1 . ■ ■. 
'. 1 

1 ■ 

• 1 

f ' 1 , ■ ; 








•./• •»• s '/••■•* /• -x 



1.0 1.5 , , 2.0 2.5 3.0 

co (a.u.) 



Figure 3: Transition probability for the spiked harmonic 
oscillator, as functions of the field frequency u> and different 
parameters £ introduced in (|79|) . We consider the transition 
from the energy level n — 2 to m = 3 to first-order perturba- 
tion theory with respect to the external laser field. The field 
amplitude is taken to be weak Eq = 0.005 a.u. and the cou- 
pling constant is chosen as A = 0.5. The pulse length r and 
the frequency u) are indicated in the figure. 

As this example is completely solvable it serves well to 
illustrate various general features. First we use it to argue 
that adding xE(t) instead of XE(t) to H in order to con- 
struct the non-Hermitian version of the Stark-LoSurdo 
Hamiltonian is unphysical. To see this we compute the 
transition amplitude in the weak intensity regime to first 
order (|5"Tj) , where instead of the amplitude ($ n |X$ m ) 
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we would have (<p n \ r\xr) 4> m ). Now for ((80)) we would 
have that {4> n \ rjxr] -1 4> m ) — (cj) n \ x (j) m ) for n ^ m, such 
that no effect would be visible in the transition amplitude 
to first order. 

Let us therefore take instead the transformation fj = 
exp(— £p 2 ). We then compute the canonical variables to 



X = x - i2£p 



and 



P = p. 



(81) 



For the corresponding non-Hermitian system H(x,p) — 
we evaluate next the transition amplitude 
for <t>2' 2 — > 4% 2 with A = 0.5, subjected to a monochro- 
matic linearly polarized electric field E(t) = Eq sin(wi) 
and depict the result in figure 3. 

As expected we obtain the main contribution for the 
transition at e^ 2 — e\ 2 = 2. The value £ = is perfectly 
reasonable and corresponds to adding XE(t) to H, with 
X given in (|80p for the reasons outlined above. However, 
for large enough values of £ we observe that the transition 
probability becomes larger than 1, which is of course in- 
consistent and unphysical. Therefore to add xE(t) to H 
is meaningless in our framework, unless x can be chosen 
to be an observable in the non-Hermitian system. 



C. The -x 4 potential 

A further interesting Hamiltonian arises when we spec- 
ify in equation {1} the parameter N — 4, which involves 
a potential which is unbounded from below. Recently 
Jones and Mateo Q established that this Hamiltonian is 
in fact isospectral to the Hermitian Hamiltonian 



H=p 2 



fx 4 



2gx for x, g € 



(82) 



This Hamiltonian is of great interest as it serves as a sim- 
plified version for the — 4 quantum field theory, which 
may for instance be used to mimic the Higgs mechanism. 
To obtain the Hamiltonian (f8"2")) from (JXJ) with N = 4 one 
needs to pass via two auxiliary Hamiltonians as follows 



H(z) Z -H H{x) -^h = rf 



(83) 



All manipulation in (|83[) arc spectrum preserving. In the 
first step the general idea @ was used to map the con- 
tour from within the wedges Wl and YVr back to the real 
axis. As discussed in section II A 1 there are many pos- 
sible parameterization, which guarantee the appropriate 
boundary condition. Unfortunately, there is no construc- 
tive method to select out the most useful contour within 
the wedges and this choice remains a matter of inspired 
guess work [8j . Here the best choice is guided by the de- 
sire to be able to construct a similarity transformation 77, 
which maps the non-Hermitian Hamiltonian H adjointly 
into a Hermitian Hamiltonian h. Hitherto, this proce- 
dure was only successful in an exact manner for the class 
of Hamiltonians in fTJ) with N = 4, in which case 77 can 
be constructed exactly either b y op erator methods Q, 
differential-equation techniques [45J or Moyal products 



[39(. Even for the next example N — 6 the same trans- 
formation used as i n p does not yield an exact similar- 
ity transformation [51| . The last step in (|53"j) in the case 
N = 4 is to transform h into the Hamiltonian H (|52")l via 
a Fourier transformation. 

Concretely, we exchange now the constant g by e in ((T|) 
with N — 4 and obtain H = —d 2 /dz 2 — ez* thereafter. 
Using now the parameterization z\(x) as defined in (JT]) 
one obtains the non- Hermitian Hamiltonian 



H* = p 2 



P , -2 



a 



{x,p 2 } 



2 ax , (84) 



The domain of H x is now the entire real axis, where a = 
16e and the coupling constant g has been introduced to 
separate off the non-Hermitian part [8|, :39j . Next we want 
to compute ij by means of Moyal products. For this we 
have to convert H first into a real valued function and 
have to substitute the anti-commutator with the Moyal 
products. Thus we have to replace {x,p 2 } by x * p 2 + 
p 2 * x — 2xp 2 . Subsequently we can use (|64]) and the 
differential equation (|57|) for the Hamiltonian ([51)1 in the 
unknown quantity Tq 2 {x,p) becomes 

= 4:gp 2 xr] 2 - 8gxar] 2 - Axad p r] 2 (85) 
- d x 7j 2 + 4:pd x Ti 2 + 2gpd p d x ii 2 - gxd 2 rj 2 . 

We can solve this by 



v 2 = e¥ ^-2 91 



(86) 



such that 77 = e e° bp. From (|2~4"|) we obtain thereafter 
the Hermitian Hamiltonian 



h x 



: V 



P 



+ a (x 2 - l) + 



9 



(p 2 -2a) 
4a 



(87) 



Let us compare how these expressions are obtained by 
means of operator identities. In principle we have to 
make a general ansatz to find q, but having already found 
77 we can simply extract it from 



(88) 



9 3 n 
q = -z-p - 2gp 
da 



and verify the corresponding expressions. From (|84|) we 
find that 



hf =p 2 - P -+a(x 2 -l) 



(89) 



Next we compute the 71-fold commutators 

c^(hf) = Aigax- ig{x,p 2 } (90) 



~(2)f/,5 



. iK)= . 
4 3 W) = o. 



(91) 
(92) 



With I — 2 we then find that the generic expression 
for the Hermitian Hamiltonian yields precisely (|87|) and 



12 



the generic expression (1601) for the non-Hermitian Hamil- 
tonian gives (|53|). 

Now the non-Hermitian system in terms of its canoni- 
cal variables 

X = x + ^- (p 2 - 2a) and P = p, (93) 
2a 

results from £P 4 (x, p) — h xi (X, P) . In addition we verify 
h x \x,p) = {H xi Y(X,P). 

In this case it suffices to choose P = p as an observable 
to make the metric unique. Note also that it is not possi- 
ble to demand X — x to be an observable as we can not 
find a function f{x) such that all functional dependence 
on p is eliminated from the term q + /(ft-p )■ 

V. CONCLUSIONS 

Given a non-Hermitian time-independent Hamiltonian 
H, we argued that the analogue of the Stark-LoSurdo 
Hamiltonian should be 

H l (t) = H + XE(t), (94) 

where X = r]^ 1 xrj is the position operator in the non- 
Hermitian system. As we have shown when we simply 
add xE(t) to H, we obtain unphysical results unless x is 
an observable in the non-Hermitian system. However, we 
also demonstrated that this is not always possible and x 
is often degraded to be a mere auxiliary variable in the 
non-Hermitian system. 

As in the time-independent scenario we saw that once 
the similarity transformation is known, one can easily 
translate all the relevant calculations into the Hermitian 
system. The situation is less straightforward when the 
transformation r\ and therefore the Hermitian system is 
not known. In that case one may take our expressions as 
benchmarks and think of various different approximation 
schemes, such as standard perturbation theory, a pertur- 
bation via the C-operator, Floquet type approximations 
for periodic potentials etc. 

From what has been said one may adopt a rather pes- 
simistic standpoint and conclude that in the end the non- 
Hermitian formulation is in most cases a mere change of 
metric of a well posed Hermitian problem. Nonetheless, 
even leaving the technical difficulty aside to establish the 
precise relation between these conceptually different for- 
mulations, it has been successfully argued that the non- 



Hermitian formulation is often more natural and simpli- 
fies computations [52l . |53| . For an atomic physicist this is 
of course a natural scenario when we compare these al- 
ternative formulations with treatments in various gauges, 
which are also just different ways to express the same 
physical quantity. It is a well established fact that differ- 
ent choices of gauges often drastically simplify problems 
in that context and allow for a more intuitive interpreta- 
tion. For instance, tunneling processes can be visualized 
and interpreted more easily in the length gauge formu- 
lation, since then one may picture the problem in terms 
of a time-dependent effective potential barrier, whereas 
all other gauges would obscure this intuitive physical in- 
terpretation. Furthermore, phenomena occurring in the 
context of high frequency fields are most intuitively un- 
derstood when viewed in a time-dependent dichotomous 
potential in the Kramers-Henneberger gauge 

Let us conclude by commenting on some of the immedi- 
ate open problems, which follow from what we discussed. 
Concerning the time-dependent treatment it would be in- 
teresting to change the current set-up by allowing r\ to 
be time-dependent. 

Having entirely focussed on the pseudo-Hermitian na- 
ture of the Hamiltonians involved, we want to conclude 
with a final comment on the role played by PT-symmetry 
in the time-dependent setting. When [PT, rj\ = the 
term XE{t) is only PT-symmetric when E{—t) = —E(t). 
This means that 'PT-symmetry depends on the explicit 
form of the laser pulse. Taking for instance a typical 
pulse for a laser field with frequency lu, amplitude Eq 
and Gaussian enveloping function /(£), that is of the 
form E(t) = E sm(ujt) f (t) , the term xE(t) would be 
PT-invariant. However, the perfectly legitimate replace- 
ment sin((x)i) — * cos(cji) in this field would break the 
PT-invariance. Recall that in this context the electric 
field is treated classically. For a discussion of PT- 
symmetry for a full quantum electrodynamic setting we 
may refer to [H], [EH|. However, for the physical appli- 
cations we dealt with in this manuscript, PT-invariance 
is not a relevant issue, since the pulse is always chosen 
such that H<S>(0) = e$(0) and H<P(t) = e$(r). The 
consequences of PT-symmetry on the eigenvalue prob- 
lem is therefore only important when considering the full 
time-independent eigenvalue problem (|32p . To investi- 
gate this full solution of (p{2"|) . the consequences on the 
non-Hermitian counterpart with its dressed states [32| 
would be extremely interesting [56| . 
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